Quasi variational solutions for first order 
quasilinear equations with gradient 

constraint 

Jose Francisco Rodrigues Lisa Santos 

Abstract 

We prove the existence of solutions for an evolution quasi- variational inequality with 
a first order quasilinear operator and a variable convex set, which is characterized by 
a constraint on the absolute value of the gradient that depends on the solution itself. 
The only required assumption on the nonlinearity of this constraint is its continuity and 
positivity. The method relies on an appropriate parabolic rcgularization and suitable a 
priori estimates. We obtain also the existence of stationary solutions, by studying the 
asymptotic behaviour in time. In the variational case, corresponding to a constraint 
independent of the solution, we also give uniqueness results. 



1 Introduction 

The initial-boundary value problem (IBVP) for scalar equations of first order 

dtu-V- *(«) = f(u) (1) 

in bounded sets Q C 1^ with a smooth boundary dQ and time t > may, in general, be 
not well posed, even for smooth flux function $ = &(x,t,u), source term / = f(x,t,u) 
and initial-boundary smooth data (here we denote by dt the partial derivative in order 
to t and by V- the divergence in x S f2). In the pioneer paper [4j, Bardos, Leroux and 
Nedelec, using the method of vanishing viscosity, extended to bounded domains the notion 
of entropy solution, obtaining their existence and uniqueness in the BV framework. Dealing 
with data merely in L°°, Otto ([21j. see also [19] or [13j ) has shown the well posedeness 
of the IBVP introducing an appropriate weak formulation. Recent results on L°° entropy 
solutions were obtained in [20] and in [12], where the delicate boundary trace question is 
analyzed and more references can be found. The need of weak formulations is justified, 
not only by the possibility of shock fronts but also, in general, by the boundary layers 
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introduced by the vanishing viscosity and the impossibility of prescribing Dirichlet data for 
solutions of dH) on the whole boundary d£l. 

Here we are interested in solutions of ([I]) with the spatial gradient globally bounded, 
which not only prevents the existence of shocks but also allows the prescription of the data 
on the whole d£l. This will be done by imposing a gradient constraint to the IBVP and 
has a motivation from critical-state problems. 

It is known that several critical-state problems arising in different physical models, such 
as in elastic-plastic deformations, sandpile growth, magnetization of type-II superconduc- 
tors or formation of networks of lakes and rivers, are formulated with convex constraints 
on first order derivatives (see, for instance, [BJ, [13] or [52] and their references). More than 
one approach has been considered leading, in particular, to several works with gradient 
constraint for the scalar function u = u(x,t), 



where G = G(x, u) > is a given threshold function. For the elastic-plastic torsion problem, 
the variational inequality formulation of the elliptic problem for the case G = 1 is well- 
known (see [31] , [9] or [25] , for instance) and has been extended to elliptic quasivariational 
inequalities in [17], when G = G(u). 

The gradient constraint ([2]), for each t, splits the domain Q into two regions 



which are, formally, separated by an unknown free boundary dl(t) flO = dA(t) n 0,. In this 
work we are interested in imposing the constraint ([2]) to solutions of first order quasilinear 
equations of the type ([]]). In fact, we look for a function u = u(x, t) that solves the equation 
(fT]) in the a priori unknown region A(i), i.e., if the gradient threshold is not attained, while, 
in the complementary region I(t), u solves the Hamilton- Jacobi equation | Vit| = G(u), with 
the prescribed boundary and initial values u = on <9S7 x (0,T) and u(x,0) = uq(x) on $7, 
respectively. 

Similarly to the elastic-plastic torsion problem, as it was shown by Brezis [8] (see also 
[25], p. 266), we may introduce a formal Lagrange multiplier A = X(x,t) associated with the 
constraint ([2]) and such that 



Vu| < G(u) 



(2) 



A(t) = {i£fi: \Vu(x,t)\ < G(x,u(x,t))} 
I(t) = {x G tt : \Vu(x,t)\ = G(x,u(x,t))} 



A > 



\{G(u) - \Vu\) =0, 



(3) 



and 



d t u - V • *(u) - V • (AVn) = f{u), 



(4) 
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where now these conditions hold in the whole space-time cylinder Q = Vt x (0, T) in some 
sense. Then, introducing for each t G (0, T) the convex set 

KgM*)) = {v e Fq 1 ^) : |Vt;| < in ft} 

and integrating by parts with — for v 6 I&G(u(t)) we observe that, for u = u(t), 
([3]) yields 

/ AVu-V(w-u)< f \(\Vu\ \Vv\ - \Vuf) < [ \(G(u) - \Vu\)\Vu\ = 
Jn Jn Jn 

and then u(t) G ^-G(u(t)) satisfies the quasivariational inequality, for a.e. t G (0,T), 

f d t u(v-u)+ [ *(«) • V(« - tt) > / /(«)(«-«), (5) 

Jn 

with the prescribed initial condition u(0) = no- 

This is a natural formulation for general first order quasilinear scalar operators under 
the constraint (J2J) and it leads to a new class of problems that were not considered in 
the classical references of quasivariational inequalities [3] or [7j- However, it contains as a 
special case the problem considered in [5] motivated by a critical-state model in supercon- 
ductivity. It was also considered in the setting of parabolic operators in |28t [29], in the 
variational inequality case, i.e., when G does not depend on the solution u, and in a qua- 
sivariational case by the authors in [2B]. We note that Prigozhin has considered a similar 
evolutionary variational inequality in superconductivity [23] and for a sandpile model in 
|24j . corresponding essentially to the particular case <& = 0, / =constant> and G = 1. 
This last model, which can also be interpreted as an evolution governed by the multival- 
ued differential associated with the constraint |Vf | < 1, has also been considered by other 
authors (see, for instance, [1] or [10] and their references). Recent results for parabolic 
problems with gradient constraint have been obtained in [2] and |15j . 

Under the natural assumptions on the smoothness of $ and / and only positivity and 
continuity on u of G, we prove the existence of a globally continuous solution u to the qua- 
sivariational inequality ([5]), with a bounded measure time derivative dtu. For the special 
case when G = G(x) > is merely bounded we show the uniqueness of solution of 
the variational inequality ([5]) in ~Kq and the L 2 integrability of dt$. Under appropriate 
assumptions on the data we also prove the asymptotic behaviour of the evolutionary solu- 
tions as t — > oo. We obtain, in particular, the existence of a solution u°° to the stationary 
quasivariational inequality ([5]), when $ is time independent, / is strictly decreasing in u 
and G is bounded in m or V • $ + / is compatible with the weak maximum principle. In 
the variational inequality case, if /°° is strictly decreasing in u, we have uniqueness of 
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the stationary solution and we obtain an estimate of the asymptotic stabilization of 
$(t) — > i? 00 as t — > oo, first in L 1 and then also in Holder spaces. 

The precise assumptions and main results are stated in Section 2. We use the method 
of vanishing viscosity and constraint penalization in two consecutive steps, by considering 
an approximating quasilinear parabolic problem with convection, extending the approach 
of [26]. The essential a priori estimates are established in Section 3 and the passage to the 
limits in Section 4. In the intermediate step, we show the existence of a solution to the new 
parabolic quasivariational inequality (jlip . generalizing a result of [26] in the case p = 2. 

The existence of unique solution of the variational inequality case is proven in Section 

5, extending the techniques used in [2S] and [22] for the parabolic case. We also give an 
indirect proof of the existence of the quasivariational stationary solution in the final Section 

6. There we prove the asymptotic behaviour as t — > oo, extending the estimate of [26] in dtu 
for the quasivariational case and adapting an argument used in the parabolic variational 
inequality setting of [TT] when the threshold does not depend on the solution . 

Finally we observe that all the results of this work also apply to linear transport equa- 
tions 

dtu + b • V-u + cu = F, 

under appropriate assumptions on the data b, c and F, but since this case presents inter- 
esting additional properties it will be treated in the forthcoming article [27] , 

2 Assumptions and main results 

Let f2 be a bounded open subset of M. N with a Lipschitz boundary dQ and, for T > 0, let 
Q = Q x (0, T). We denote by x = (x±, . . . ,xn) an arbitrary point of 0, V = (d Xl , . . . , d XN ) 
denotes the spatial gradient and dt the partial derivative in order to the variable t G (0, T). 
We represent vectorial functions using a bold symbol, F = (iq, . . . , Fn). Its divergence is 
then V ■ F = d Xi Fi = divF. 

Assume $ = &(x,t,u) :QxE-> R N , f = f(x,t,u) : Q x R — > R are such that, for 
any R > 0, 

* G W 2 '°° (Q x (-R, R)), f e W 1 ' 00 (Q x (-R, R)) , (6) 

and V • $ and / are functions with at most linear growth in the variable u, uniformly in 
(x,t), i.e., there exist positive constants c\ and C2 such that, for a.e. (x,t) E Q and all 

|(V • *)(x,i,-u) + f(x,t,u)\ < ci\u\ + c 2 . (7) 
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For the gradient constraint G = G(x, u) : O x R -> R we assume 

Ge C(R;L°°(ft)), G{x,u)>\>0 a.e. x G ft, Vu G R, (8) 

in particular G is uniformly continuous in u and bounded in x in each compact of R. 
Given a positive function (p G L°°(ft), we define 

= {v G : |Vf| < ip a.e. in ft}. 

Theorem 2.1 Under the assumptions (0J)-([2P, /or eac/i 

w G K G(uo) n C (ft), (9) 

i/iere exists at least a quasivariational solution 

u G L°°(0,T; Wo' 00 (ft)) nC(Q), S^tt G L°°(0, T; M(ft)) , 

to t/ie problem 

u(t) G IK G («(i)), ^(0) = u , and for a.e. t G (0,T), 

- + /" *(u(t)) • V(u - u(t)) (1Q) 



> / f(u(t))( V -u(t)), VveK G{u{t) y 

Jn 

Here M(ft) denotes the space of bounded measures in ft, the first integral in (jlOp is 
interpreted in the duality between the spaces M(ft) and Co (ft), which denotes the space of 
continuous functions on ft which vanish on <9ft. The other integrals are interpreted in the 
usual Lebesgue sense in ft. 

We shall prove this theorem using the vanishing viscosity approach. In particular, we 
shall obtain the existence of solution to the following parabolic quasivariational inequality 
with nonlinear convection and nonlinear reaction terms, extending a result of [26], 

( u s (t) G K G{uS{t)) , u s (0) = u s and for a.e. t G (0,T), 

d t u\t)(v-u s (t)) + 6 [ Vu s (t)-V(v-u s (t))+ [ $(u s (t))-V(v-u s (t)) 



(11) 

> / f(u d (t))(v-u\t)), V«eK G( „ i(t)) . 
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Proposition 2.2 Under the assumptions |7|)-([3)) and if u 5 Q G ^g(m ) ^ Co(^) is such that 
A-Uq G Af(f2), £/iere exists a function 

u 5 g l°°(o,T; w 1,oo (^)) nC(Q), G L°°(0,r;M(fi)), 

which is solution to the parabolic quasivariational inequality Ml}) . 

When the function G does not depend on u the conclusions are stronger and we can 
obtain the uniqueness of solution in the corresponding problem, which is then a variational 
inequality. 

Theorem 2.3 Under the assumptions {6]j, ([?]) and 

< A < G(x) < A for a.e. x £ Q, (8') 
n G Kq, 

there exists a unique function 

# G L°°(0, T; Wq 1 ' 00 ^)) Pi C(Q), d t $£ L°°(0,T;M(n)) nL 2 {Q), 

such that •& is the solution of the first order variational inequality 
( G Kg, 0(0) = ii and /or a.e. t G (0,T), 

(12) 



a t t?(t)(v-i?(t)) + / *(#(*))• v(u -#(*)) 

f2 Jf7 



> / /(0(i))(v-0(t)), Vv£K G . 



We prove the stabilization in time (for subsequences t n — > oo) to steady-state solutions, 
both in the quasivariational and in the variational cases. Given functions f°° = f°°{x,u) : 
UxR — ► R and * = u) : x R — ► R^, we consider the stationary quasivariational 
inequality: to find 

u°° G K G ( u oo) 

f f (13) 

Theorem 2.4 Assume (EJj, 3> is iime independent and f is strictly decreasing in 

u, i.e., 

d t ® = o, d u f < -« < o (14) 
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and 

< A < G(x, u) < A for a.e. x G Q, for alluGR (15) 
or there exists M > smc/i £/ia£ ; /or aZZ i? > M, 

(V • *)(sc,12) + f(x,t,R) < 0, (V • *)(x,-i?) + f(x,t,-R) > 0. (16) 

Setting £ji(t) = / sup |3t/(x, t, u)| cZrr, we suppose, in addition, that for R>Rq, 
Jn \u\<R 

ft+l rt+l 

sup / tn{T)dT < C R and / CR(r)dr — > 0, (17) 

0<t<oo Jt 



feeing C^j a positive constant and 

f(x,t,u) — > f°°(x,u), for all \u\ < R and a.e. x G ft. (18) 



t— >oo 



Then problem h!3\) has a solution which is the weak-* limit in Wq'°°{VI) and strong 
limit in C 0,a (ft), < a < 1, for t n — > oo, of a sequence {u(t n )} n , being u a global solution 
of problem MCfy . 



These asymptotic results still hold when G is independent of u. However, in this 
case, when we can guarantee the uniqueness of the stationary solution for the variational 
inequality, we may conclude the convergence of the whole sequence, i.e., 

•&{t) in W^°°{Q) weak- *, 

and, by compactness, also strongly in spaces of Holder continuous functions. 

Proposition 2.5 Suppose that 3>, f°° satisfy |7^P , (EJ), (7p, G = G(x) satisfies [8]) . 
Then the following stationary variational inequality 

( $(tf°°) • V(u-i?°°) > / Z 00 ^ 00 )^-^ 00 ), WgK g , 
Jn Jn 

has a unique solution 

With additional assumptions we may even estimate the order of convergence of $(t) to 



7 



Theorem 2.6 Under the assumptions of Proposition [275\ and 

rt+1 r 

1 M (t) = / / sup \f(x,T,u) - f°°(x,u)\ dxdr — > 0, 
Jt Jq \ u \<m 

for M > Ma, i/i? denotes the solution of the variational inequality {7j|] andil) 00 the solution 
of the variational inequality [T3\) then 

#(t) — > in C°' Q (fi) /or < a < 1. 

t— >oo 

In addition, if rj M (t) = 0(e~ 7 *), 7 > 0, lo/ien i — > 00, f/ien, for any < a < 1, 

\\m - 0°°|lc<>.«(fi) = 0(e"^), as t -> 00, 
where 7' = (^77)^ wi/i ^ = min{//, 7}. 



3 The a priori estimates for the approximating problem 

For < e < 1, < 5 < 1, we consider the approximating quasilinear parabolic problem for 
u eS (denoted, for simplicity, also by w), 

d t w-V- [5k £ (\Vw\ 2 - G 2 E (w))Vw + = /(to), in Q, 

( 2 °) 

io(0) = in H) = 0onffix (0, T), 

where fc £ is a smooth function such that k E (s) = 1 if s < and & e (s) = if s > e, 
G e = G*p e >A>0, being /3 e a mollifier in (x, it) and uff an approximation of uq, 
belonging to such that |Vuq 5 | < |Vtto| and || Au^H^im) < Co/5, where Cq is a 

constant independent of e and 5. 

General parabolic theory for quasilinear non-degenerate equations ( \16\ [T8] ) yields the 
existence of a unique solution u eS 6 C 2 (Q) n C(Q) to (f20|). 

We obtain several a priori estimates for the solution u eS of the approximating problem, 
independent of e and 5. 

Proposition 3.1 Under the assumptions (Ej)-(0|), the solution u eS of problem \20\) is such 
that, 

1. for a.e. (x, t) G Q, 

\u e5 (x, t)\ < M = inf e AT J max|u | + 1, f — ^— V L (21) 
A>&! o \A — 61/ 
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where b\ = c\ + \ and 62 = c i an d c 2 positive constants defined in §7ty. In 

particular, M is independent of e and 5; 

2. there exist positive constants C\, C2 and K , such that 

!MU(o,T;L 1( n)) < e^ T {C 1+ K\\Vu^\\ LHQ) y, (22) 

3. there exists a positive constant C such that 



\k £ (\Vu £& \ 2 ~ G^{u £5 ))\\ Loo{ ^ T . L1(n)) < pC; 



4- for all 2 < q < 00, there exists a positive constant D q such that 

where D q depends on q and is independent of e and 5. 
Proof of Proposition 13. lM . 

By Theorem 2.9 of [18] . page 23, if a function w is a classical solution of the equation 



d t w - aij(x,t,w,Vw)d XiXj w + a(x,t,w,Vw) = 0, (24) 
with a,ij and a bounded and satisfying the following two conditions, 

<2jj(x, t, w, 0)£i£j > 0, w a(x, t, w, 0) > —b\w 2 — 62, 
for all £ = (£1, . . . , £at) £ M. N , with 61 and 62 nonnegative constants, then 

AT 



max t/)(i, t) < inf e < max \w 



Q 1 1 a>6i I nx{o}uanx(o,T) — &i 

The equation in (f20|) is of the type of ([2"4l . with 



ay(a;, Vto) = 25^(|Vu;| 2 - G?H)S Xl «; d^w + <5£; e (|W;| 2 - G*M)<Ji 



o(x,t,«),Vffl) = 2<5^(|Vu?| 2 - G e (w))G e (w) [d u G £ (w)\Vw\ 2 + VG £ (w) ■ Vw] 

- (V • *)(x,i,w) - (d u $)(x,t,w) ■ Vw - f(x,t,w) 

and 

w a(x, t, w, 0) = —w ((V • &)(x, t, w) + f(x, t, w)) . 



By (J7J), |(V • &)(x, t, w) + f(x,t,w)\ < 2 c\\w\ + 2 C2 and so wa(x, t, w, 0) > —biw 2 — b%, 
for 6i = 2ci + 1 and 62 = Cn. We also have \w(x,t)\ < max|tio| + 1 for a.e. (x,t) £ 
ft x {0} U dO, x (0, T) and the conclusion follows. □ 

Proof of Proposition 13. l\ -2. 

Setting v = dtw and differentiating the first equation of ([20]) in order to t, we get 



d t v 



6k' £ (\Vw\ 2 - G 2 £ {w)) (2 8 Xi w d Xj wd Xj v - 2 G E {w){d u G E ){w) v d Xi w) 



+5k £ {\Vw\ 2 - G 2 (w))d Xi v + (dt$i)(w) + (d u $i)H<> = (dtf)(w) + (25) 

We define the signal function by sgn°(0) = 0, sgn°(r) = 1 if r > and sgn°(r) = — 1 
if t < and we consider the approximation by a sequence of C 1 increasing functions, 
s c : R ->■ R such that 

C 

< s' c ( r ) < -T for a11 T G M > s d T ) = s g n °( r ) for a11 t £{0}U R\] - C, C[- (26) 



Setting S^(t) = / s^(a)da, we have 
■/ 



lim ts'At) = 0, lim SAt) = Irl, Vr G R. (27) 

As f = on dQx]0, t[ we also have s^(v) = on dQx]0, t[. If we multiply ([25]) by s^(v) 
and integrate over = $7 x (0, t), we obtain 



d t vs c (v) + 2S / ^(IV^P-G^H)^^^^^^^^^^) 
■/Qi 

-26 [ K(\Vw\ 2 -G 2 £ {w))G £ {w){d u G £ ){w)vd Xl wd Xl vs' (: {v) (28) 



+ 5/ k £ {\Vw\ 2 - G 2 £ (w)) \Vv\ 2 s'Av) - (a Xi dt$i)(w)sc(v) 
Qt JQt 



(dtf){w)s c (v)+ / (3J)H« ( W. 
JQ< ■/ Qt 

We note that 

0) = d t w(x, 0) = SAu e s + (V • *)(*, 0, uf) + (d u &)(x, 0, u £ s ) ■ Vuf + /(x, 0, uf). 
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Let us analise each term of the left-hand side of the equality ([28 



Qt 



Mv)= / [S c (v(t)) - S c (v(0))] = / S c (d t w(t))- / S c (d t w(0)) 
at in Jn Jn 



and 



lim / dtv s^(v) = / |(9tw(t)| 
^ JQ t Jn 

- ! + (V • *) (x^uf) + (d u &)(x,0,u e o) -V^ 5 + f(x,0,u E s )\. 

Jn 

On the other hand, 



25 / k' £ (\Vw\ 2 -G 2 £ {w))d Xi wd Xj wd Xj vd Xi vsUv) 
JQt 



25 / k' £ (\Vw\ 2 -G 2 (w)) (VwVv) 2 s' ( {v) >0, 



and 



5 k £ (\Vw\ 2 -G 2 (w))\Vv\ 2 s[(v) >0, 



< if, 



< if ||Vw|| L i W) 



where, for M defined in (f2"T|) . if = ||^||vk 2 >°°(Qx(-m,m)) is independent of e and (5. By ([2 
and Lebesgue Theorem 



and 



-2 5 / femV^p-G^HjGeH^GeJH^wS^uus^t;) — ► 0, when C ^ 
JQt 

/ (d u <&)(«;) • VvvsUv) — > 0, when £ ->■ 0. 
Gathering all the information above, 

f \v(t)\< [ \5Aul 5 + (V ■^)(x,0,uf l s ) + (d u ^)(x,0,uf ) s )-Vul 5 + f(x,0,ul 5 )\ 
Jn Jn 

+ K\\Vw\\ Ll(Q) + J \dtf(w)\ + J \d u f(w)v\, (29) 



and we obtain 



j \v(t)\ < d + if ||V«;|| il(Q) + C 2 




v . 



o Jn 
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Here, the constant C\ is dependent on ||tto||flJ(n)> ^II^^IU 1 ^) (bounded by Co), 
ll^llw^^CQxC-Af.M))' \\f\\w^(Qx{-M,M)Y the constant C 2 depends on \\f\\ w ^(Qx(-M,M) 
and K depends only on ||3?||w 2 .°°(Qx(-M,.M'))- The conclusion follows immediately by ap- 
plying the Gronwall inequality. □ 

Remark 3.2 When $ is independent oft, we may take K = in the last theorem. 



Proof of Proposition 13. 11 -3. 

We multiply the equation (|20p by w and integrate over S7, obtaining, for a.e. t E (0, T), 

5 [ k £ (\Vw(t)\ 2 - G 2 £ (w(t)))\Vw(t)\ 2 = [ (f( w (t))-d t w(t))w(t)- [ 3>(w(t))-Vw(t). 
Jn Jn Jn 

So we get for w = w(t), 

5 [ \Viv\ 2 <5 [ k £ (\Vw\ 2 -G 2 (w))\Vw\ 2 
Jn Jn 

1 II 2 8 [ I ,2 

<M\\f-d t w\\ Loo{0T , Lim + — ||*H|| L » ( o iT . ia( n)) + 2 ' \ Vw \ > 



from which we conclude 



/ \Vw(t)\ 2 < [ k e (\Vw(t)\ 2 -G 2 (w(t)))\Vw(t)\ 2 < for a.e. i E (0, T). (30) 

As k £ (s)s > for s > 0, fc e (s) = 1 if s < and G e > A, we obtain the following estimate, 
for a.e. t G (0, T), 

A 2 / M|V^t)| 2 -G E 2 Mt)))< / A; £ (|Vu;(i)| 2 -G 2 Wt)))G 2 Ki)) 

Jn 



k £ (\Vw(t)\ 2 - G 2 (w(t))) (G 2 (w(t)) - \Vw(t)\ 2 ) 

+ / k £ (\Vw(t)\ 2 -G 2 (w(t)))\Vw(t)\ 2 
Jn 

< I G 2 (w(t)) + ^C< {\\G\\ 2 LX + 1) \Q\ + 1(7. 



Proof of Proposition 1331 -^ . 
Calling 

A £ (t) = {xen-. \Vw(t)\ 2 > G 2 (w(t)) + e}. 



□ 



12 



we decompose Q = A £ (t) U (fi \ A £ {t)). Then 

\Vw(t)\ q < [ (G 2 £ (w(t))) + e)" < (\\G\\ 2 LOO + l)" (31) 
o\A s (t) Jn 

Recalling ([30]) and, as k E (s) = e~ > -f^, for s > e, we have, for g even integer greater 
or equal to 4, 

(\Vw(t)\ 2 - G 2 (w(t)))^ . ,. |2 



\Vw(t)\ 



k £ (\Vw(t)\ 2 - G 2 (w(t))) \Vw(t)\ 2 < (32) 



M*) 



Then, using (|32j) we find 



9-2 



\Vw(t)\ q = / [f|Vu;(t)| 2 -G 2 (^(t))) +G 2 (w;(t))l 2 |Vu;(t)| 
AM) JAM) LV ' J 



As(t) 



9-2 



< / 2^-M(|V^(t)| 2 -G 2 Kt))) 2 +Gr a (w(t)))|Vu>(t) 



< 2 ^ e ^(t2)! / fc e (|V^(t)| 2 - |V W (t)| 2 (33) 

JAM) 



r(f) 

+ 2^ / G|- 2 (io(t)) |Vio(t)| 2 

•M £ (t) 

9 ~ 4 II || <7— 2 || 1 1 2 
-£2 C 9 + 2 2 II^IIl 00 ll Vl,j llL°°(0,T;L2(n))- 

Using again (i30"j) . (f3Tj) and (j33j) . we obtain that 

II V W llL°°(0,T;L9(f2)) - ^2 1 



□ 



4 Existence of solutions by letting £ — > and £ — >• 

In this section we prove Proposition l2.2l and Theorem l2.ll We start with a simple lemma: let 
a sequence {h £ } of functions belong to L°°(r2) and converge to h, where h and h e are greater 
or equal to A > and define a 6 = \\h — /i £ ||l°o(q). Then, for any v G K^, v e = x _^ ae v G IK/^ 
and u e — > i> in Wq' 00 ^). This strong approximation of a given function belonging to 
by a sequence of functions belonging to K^e will be a key tool in the proof of existence of 
solution to the problems (fTTI) and (fTOl) . 
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Proof of Proposition 12.21 

By the estimates (|22[) and ()23|) . we have the uniform (independent of e) boundedness of 
{d t u eS } £ in L^^TjL 1 ^)) and of {Vu e5 } £ in L°°(0, T; L q (Q)). This implies the conver- 
gence, at least for a subsequence, 

d t u £& d t u s in L°°(0,T; M(Q)) weak- * when e -»• 

Vu rf Vu 5 in L°°(0, T; L 9 (0)) weak- * when e -> 0. 

By the compactness of the imbedding W Q ' q (^l) C C(Cl), for q > N, and a theorem on 
compactness of functions with values in an intermediate Banach space (see [30], page 84) 
{u eS } is relatively compact in C([0,T];C(fl)). So, there exists a subsequence such that 

u £& — >u s in C{Q), whene^O 

and, consequently, G £ (u £& ) — > G(u & ) in L°°(Q) and u s (0) = Un in O. 

We are going to prove that u 5 (t) G K G / u 5/ t ^, for a.e. i G (0, T). Set ^4 e = {(x,t) G Q : 
| V« £l5 (x, t)| > G £ (u eS ) + v^}- Then, by the definitions of /s e and the uniform boundedness 
(in e) of {£; £ (|Vu £<5 | 2 - G 2 (u e5 ))} £ in L X {Q), 

f f k £ (\Vu £S \ 2 - G 2 (u £S )) i 
\A £ \ = 1 < / 1 x eV ^ < Ce 

and this implies that \A F \ — > 0. So 

[ (\Vu 5 \ -G(u 5 )) + = [ liminf (\Vu e8 \ - G(u eS ) - ^) + 
JQ JQ £ 

< liminf f (\Vu £5 \ - G(u £& ) - y/e) X 

Jq 

< lim |||Vu £<5 | -G{u £6 )\\ T2ln , \AA* =0 

and \Vu s \ < G(u s ) for a.e. in Q, i.e., G ^G(u s (t)i f° r a - e - * e (0>^)- 

As we stated at the beginning of this section, given t/ G L°°(0, T; C(O)) such that 
G K G{u s {t)) for a.e. t G (0,T), we may find v eS G L°°(0, T; C(O)) such that v eS (t) G 

KG,( u rf( t )) and v eS — > v 5 in L°°(0, T; W 1,oo (fi)), since G £ (n £5 (t)) — » in L°°(0). 

Multiplying the first equation of (|20p by and integrating over f2 x (s,i), with 

< s < t < T, we obtain 



I f d t u £& (v eS -u 5 )+5 I f k £ (\Vu £5 \ 2 - G 2 {u £5 )) Vu £S • V(v eS - u 5 ) 

Js J Q, Js J Vt 

+ f f ®{u eS ) -V(v £S -u 5 ) = f I f{u £& ) (v eS -u 5 ). 
J s J q, J s Jq 
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By the monotonicity of k e and since v £ (r) G K Gs / u e«/ r ^ for a.e. r G (0, T), we have 

f Vv eS - V(v £S -u 5 ) > I I k £ (\Vu e5 \ 2 - G 2 £ (u £5 ))Vu £S ■ V(v e5 -u 5 ), 
J s Jn JsJCl 

and the limit when e — > yields 

f [ d t u 5 (v 8 -u s )+5 f f Vv s ■ V(v 5 - u 5 ) 

+ f I *(u s ) • (v s - u s ) > f I f(u 5 ) - u 5 ). 

Js Jn Js J ll 

Using Minty's Lemma and since s and t are arbitrary, we obtain as in |26| . 

I d t u s {t) (v s (t)-u 5 {t)) +5 [ Vu 5 {t)-V(v s (t)-u s {t)) 
Jn Jn 

+ [ *(u s (t))-(v 5 (t)-u s (t))> [ f(u s (t)) (v s (t)-u S (t)), 
Jn Jn 

for a.e. t G (0, T), which concludes the proof. □ 
Proof of Theorem 12.11 

The solution u s of problem (HU satisfies |Vu 5 | < G(u s ) < G* a.e. in Q and therefore, by 
assumption ([8]) and estimate ([2T]h has a gradient uniformly bounded in 5, i.e., 

\\ uS \\L^(0,T;W^ x (n)) - C ( 34 ) 

where C is a constant independent of 6. 

Letting e — > in the estimate ([22]) , and using (fM|) , we obtain also 

IML~ ( o, W n)) ^ + K \M LHQ) ) < Cs- 

So, there exists a function u such that, when 5 — > 0, 

d t u s d t u in L°°(0, T; M(fi)) weak- *, Vu in L°°{Q) weak- *, 

and, again by the same compactness argument used in the previous proof, we also have the 
strong convergence 

u 5 — >u in C(Q), 
S— >0 

which, in particular, yields u(0) = uq in J7. 
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Then, as before, for a fixed t G (0,T), given v G L°°(0, T; C(fi)) such that u(t) G 
K G*(n(t)), w e can find v 5 (t) G L°°(0, T; C(0)) such that G K G(u «( t) ) and t;" 5 — > w in 

L°°(0, T; W 1,oo (fi)) and we can take the limit 5 -> in 

r / cV -u 5 )+8 f I Vu s • V(u* - u 5 ) +/"*/" • V(u* - u 5 ) 

is 7ri J s Jq J s Jo, 

> f f /(«') (v S -v 6 )- 
Js Jn 

As in previous step, by the arbitrariness of < s < t < T, we conclude, for a.e. 
t G (0, T), 



dtuit) (v(t)-u(t))+ / *( U (t)).(u(t)- U (t)) > / /(«(*)) (u(t)-u(t)), a.e. t G (0,T). 

Jf2 J!) 

Finally, given any measurable set u C Q, we have 

/ |Vu| < liminf / |Vu 5 | < lim f G(u s ) = f G(u) 

and consequently, for a.e. t G (0, T) we have 

\Vu(t)\ < G(u(t)) a.e. in O, 
which means that it(t) G Kg( tt w) for f G (0, T), concluding the proof that u solves (fTQ|) . □ 



5 The variational inequality 

When the gradient constraint G is a function independent of u, the first order quasivaria- 
tional inequality becomes a variational inequality. As a consequence, the solution is unique 
and more regular. 

Proof of Theorem \2.3[ 

We only need to prove that dt'd G L 2 (Q) and the uniqueness of solution. Multiply the first 
equation of problem (j20|) by dtw and integrate in Q. Then 

f \d t w\ 2 + 5 [ k £ (\Vw\ 2 -G 2 £ )VwVd t w+ [ ®(w)-Vd t w= [ f(w)d t w. (35) 
Jq Jq Jq Jq 

We remark that 

/ &(w) ■ Vd t w = — I (V • <&) (w) d t w — / d u <b(w) ■ Vw d t w 
Jq Jq Jq 
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and we obtain 



I (f(w) + (V-&)(w) + d u $(w)-Vw)d t w <\ I \d t w\ 2 + C'\\Vw\\ L 2 (Q) , 
JQ J 1 JQ 

where C is a constant depending only on M and on the assumptions on <1> and /. 

Let K e (s) = / k £ {r)dT. We have K £ (\Vw(0)\ 2 - G 2 ) = \Vu s \ 2 - G £ and also 
Jo 

-K £ (\Vw{t)\ 2 - G 2 ) < G 2 £ a.e. in Q, 

since K £ (\Vw(t)\ 2 - G 2 £ ) = \Vw(t)\ 2 - G 2 £ > -G 2 £ if \Vw(x,t)\ 2 < G 2 {x) and otherwise 
K £ {\Vw(t)\ 2 -G 2 ) >0. 

^From (|35p we then obtain, 



[ \d t u e5 \ 2 < 2C'\\Vu e5 \\ L 2 {Q) + 6C". (36) 
JQ 

Letting first e — > in (|36p . we have 

f \d t u 5 f < 2C"||Vn 5 || i2(Q) + 5C" < 2C'\\G\\ L 2 {Q) + 5C" . 
JQ 

Finally, letting 5 tend to zero, we conclude that dt'd G L 2 (Q). 

To prove the uniqueness of solution we consider a sequence s<; : R — )• R of C 1 increasing 
odd functions, approximating pointwise the function sgn°, defined as in (|26|) . If $i and i? 2 
are two solutions of the variational inequality (|12j) we define ui = i?i + C 2s c($2 — for C 
small. As 

|V«i| = |Wi + C 2 4(^2 -i?i)V(i? 2 - #i)| 

= j (1 - C 2 <(^ 2 - Vi?i + C 2 <(^2 - #i)W 2 | 
< (1 - C 2 4(^2 - G + C 2 <(tf 2 — G = G, 

since 1 — C 2 s^(^2 — #1) > 0, we have v\ € Kq. 

Substituting v\ in the variational inequality (|12p satisfied by #1, we get 



> / /(^i(t)) SC (^(t)-^i(t)). 



<2 
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Using now u 2 = ?9 2 + C 2s c(^i — ^ 2 ) as ^ es ^ function in the problem for 192) we obtain 
the following inequality for u = i?i — $ 2 , 



< / (/(*(*)) -/Oh(t))H(«(*))- (37) 
Jn 

By assumption (jSJ) there exists > such that \$>(x,u) — $f(x,v) < L$|n — u|. So, 

by m 



s' c (u(t))(*(0i(t)) - • Vu(t) < L # / |V«(t)| sUu(t)) \u(t)\ — »• 0, 

applying the dominated convergence theorem and lim rs ( -(r) = 0. For some positive con 



r->0 



lim / S c (u(t)) = / \u(t)\ <L f I I \ 
^°Jn Jn Jo Jn 

The uniqueness follows by the Gronwall inequality. 



stant Lf > 0, we also have |/(t>i) — f{v2]\ < Lf\v\ — u 2 |, and so 

(/(tfi(t))-/(tf 2 (*))K(^)) / Ht)\. 

n Jn 
Then, integrating ([3T|) between and i and letting Q — > 0, we obtain, recalling ([27)) and 
fi(0) = 



□ 



Proof of Proposition I2.51 

The existence of solution to (|19p is a consequence of Theorem I2.4I since the assumption 
(fT5)) reduces to ((85 when G does not depend on it. 

The uniqueness can be shown as in the evolutive case: suppose that $1 and #2 are two 
solutions of (|19j) ; remarking that if denotes a C 1 approximation of the sgn° function as 
in (|26p . we have, in L p (r2), 1 < p < 00, 



s<(01-0 2 ) — ► sgn°(#i - tf 2 ) 



1 in > tf 2 }, 
in = tf 2 }, 
-1 in {t?! < tf 2 }. 



Using u = tfi - C 2 sc(^i - #2) m (HU) for t?°° = #i and u = i? 2 + C 2 «c(^i - #2) in CE 
for $°° = # 2 , we find, recalling that 3> is Lipschitz in u, 

(/"W - /°°(tfi))s C 0?i - 2 ) < / - tf 2 ) - 2 | G. 
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Taking the limit ( -> 0, we obtain 

(r(#2) - rw) - / (/°°(^) - rw) < o 

{#1># 2 )} V J J{Vl<t>2)} V J 

and, since /°° is strictly decreasing, we conclude that 

\{&x >#2)}| = |R <^ 2 )}| =0 
and so $i = i?2 a.e. in Q. □. 

6 The asymptotic behavior in time 



Proof of Theorem \2.4\ 

We consider again the approximating problem (j'20|) defined in Q x (0, oo) in order to show 
that dtu eS (t) vanishes uniformly in e and 8 in M(f2) as i -> oo. 

We multiply ([25]) by e^s^ with defined in ([26]) and satisfying (j2"T|) and v = 

dtu eS (t), using integration by parts in x (s,t), < s < t: 

^ T dtv(T)s c {v{T))dr = e^ [ S(;(v(t))-e^ [ S ( {v{s)) - f* [ f ief MT S c (v(t))o!t. 

Jn Jn Js Jn 



e 

s Jn 



Similarly to ()29|) . where now K = since is independent of t by (|14p . after taking 
C — > we obtain 



/ <e^ s / |d t u £(5 (s)| + /"* / ^ T |a t n e<5 (r)|dr 

Jfl Js Jfl 



+ 11 e^\d t f{u e5 (r))\dr+ I [ e» T d u f{u eS (r))\d t u eS (t)\oIt. (38) 

Js Jn Js Jn 

By the assumption (fbi|) . 9 U / < — fi < and setting R such that, 

|u e *(M)| < R V*>0, a?€fi, (39) 



uniformly in e and 5, we deduce from (J38J) that 

e^*i/(t) < e" s u(s) + ^ e^ r 6*(r) dr, t > s > 0, (40) 

where we set u(t) = / \dtu eS (t)\ and = / sup \dtf(u)\. 

Jn Jn\u\<R 
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^From (|40p and with an elementary estimate of its last term, we obtain the following 
inequality for any t > 0, a > 0, 

ps+l 

v{t + a)< e~^v{a) + sup / &(t) dr 

CT<S<t+<T J s 

with C M = (1 — e _/1 ) _1 + 1. Consequently, from ([TTJ) , we conclude 

u{t + a)<e-^{u^) + C^C R )+C^ sup [ S+ \ R (r)dT, 

cr<s<t+a J s 



and the right hand side vanishes, as t — > oo, uniformly in e and 5, provided (|39|) holds. 
But the L°° estimate of n globally in time holds with R = R(A,diam(Q)) in the case of 
assumption (|15p or with i? = max{M, ||wo||£°°(n)}, where M is given by the assumption 
(|16p in the second case. Notice that, when (|16p is satisfied, by the weak maximum principle 
(see, for instance Theorem 1.6 of [H]), R is a supersolution and — i? a subsolution of (|20p. 

Therefore, it is clear that, by compactness, we may take a subsequence t n — > oo and a 
function n°° E W 1,oo (ft) such that (0 < a < 1) 

— »• in W n 1,00 (0) weak- * and in C°' a (n) strong, (41) 

n 

d t u(t n ) — > in M(fi). (42) 

n 

For an arbitrary measurable subset w C !1, we have 

I |Vu°°| <liminf f \Vu{t n )\ < lim / G{u{t n )) = [ G(u°°), 

Jul n Ju) n Jul Jul 

we obtain that 

|V«°°(x)| < G{x,u°°(x)) a.e. in x e Q, 

and so u°° E K G ( M oo). 

Consider the quasivariational inequality (|10p at t = t n 



d t u(t n ) (v - u(t n )) + / *(«(*„)) ■ V(« - u(t n )) 

f(t n ,u(t n ))(v-u(t n )), V«GlG( u (i n )). (43) 



> 

Given er e (0, 1) and w°° an arbitrary function of K G ( UOO ) in L°°(Q), as G(u(t n )) 
G(u°°) uniformly, there exists no £ N such that, for n > no we have a.e. 

G(n°°) < _L_G(u(tn))- 
(1 - crj 
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Defining 

if{x) = (1 - a)w°°(x), 

we have 

|V« CT | = (1 - a) \Vw°°\ < (1-a) G(u°°) < G(u(t n )). 

which means that v u G K G ^ tn ^, for all n > no- So, taking v = v a in P3|) and letting 
i n — > oo we obtain 

/ *(u°°) • V((l -a)w°° -u°°) > [ / 00 (u 00 )((1-ct)w; 00 -u 00 ), 
Jn Jn 

since by the assumption (|18p and the convergence (|41[) we have 

f(t n , u(t n )) — ► Z 00 ^ 00 ) in LP(fi), 1 < p < oo. 

n 

Finally, letting a — > 0, solves problem as ui°° is an arbitrary function of 



Proof of Proposition \2.6[ 

Let be the solution of problem f)12|) and the solution of problem (|19p . 

Recall that u = i? + C 2 s c (t?°° - 1?), < £ < 1, belongs to K G . So, using v as test function 
in (|12p and calling $ = — i?, we have 



c 2 / dtmsd#) + c / <(#)*(#)• w > c 2 / 

Jn Jf2 Jn 

Using io = + C 2 sc("^) e as test function in (fT9|) . we get 

-C 2 / s' c (^)*(i?°°) • W > -C 2 f / 00 (7? 00 )s c (^). 

Jn 

Summing the above inequalities we get 

Jn Jn Jn 

After multiplication by e^ T and integration in order to r, between cr and t, since 
* _d ' ' ' rt 

dr 



e» T — / dt# 8 C (#) dr = e!* / 5 c (i?(t)) - / S c (#(a)) - // / / e^ r 5 c (i9) dr, 
n Jn Jn Jo- Jn 



we have 



n Jn Jo- Jn 



t 

Mr-t) 



[ s' ? (??)(*(?? 00 )-*(^)) -Wdr + / ^ T "*) / |/(i9°°)-/(T?)|dr. 
Jn Jo Jn 
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Using the strict decreasing property (14), we find 



jit 



a Jn 



>(S c (tf)-tfs c (tf))dT 



< 



0. 



On the other hand 



[ e r "' f b£(0)(*(0°°) -*(#)) dr< / L*|i?|4(i9)|Vi9|dr, 
Jit io 7n 

where is a Lipschitz constant, also converges to zero when £ — > 0, since lim r S/-(r) = 0. 

Therefore, we obtain 

/ |5(t + cr)| <e-"* / p(<r)|+ e^-*-^ [ \W°°)-f(0)\dT. 
Jn Jn Jo Jn 

Arguing as in [IT], page 283, we have 

||^(t) -^°°|[ il(n) < e-^||^(ff) - *°°|| il(n) +C sup ^(r)!, 

v ; v ' a<T<t 

and so 0(t) — ► in L x (fi), since ??..(£) — ► 0. As ||V(i?(t) - $°°)\\ L °° (n) < 2A then 

t— >oo *— ^ — 1 



t— >oo 



— >• in C°' a (fi), for < a < 1, by compactness. 



If, in addition, r/ M (t) = 0(e 7i ), 7 > 0, when t — > 00, then 



O(e^), 



where v = min{/i, 7} and, by a Gagliardo-Niremberg interpolation inequality, we conclude 
as in [25], page 200, for any < a < 1, 



IC°> a (f2) 



0(e-(^)^, 



as t — >• 00. 



□ 
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